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ABSTRACT 

Time dependent numerical simulations of the Ka complex of Fc XXV arc carried out as 
a function of temperature/density/radiation field variations in high-temperature astro- 
physical and laboratory plasmas. In addition to several well known features, the tran- 
sient and steady state spectra reveal the effects due to (a) time-dependent thermal and 
non-thermal radiation fields, (b) photo and coUisional excitation and ionization, and 
(c) high densities, on the 'quartet' of principal w,x,y,z lines, and dielectronic satellites. 
The highly detailed models show precisely how, assuming a temporal-temperature cor- 
relation, the X-ray intensity varies between 6.6 - 6.7 keV and undergoes a 'spectral 
inversion' in the w and z line intensities, characterizing an ionization or a recombina- 
tion dominated plasma. The dielectronic satellite intensities are the most temperature 
dependent features, but insensitive to density variations, and significantly contribute 
to the Ka complex for T < 6.7 keV leading to asymmetric profiles. The 6.7 keV Ka 
complex should be a potential diagnostic of X-ray flares in AGN, afterglows in gamma- 
ray bursts, and other non-equilibrium sources with the high- resolution measurements 
possible from the upcoming mission Astro-E2. It is also shown that high electron den- 
sities attenuate the line intensities in simulations relevant to laboratory plasmas, such 
as in inertial confinement fusion, laser, or magnetic Z-pinch devices. 

Key words: Atomic Processes - Atomic Data - X-ray: K-shell - line:profiles - 
galaxies iSeyfert - X-rays:galaxies 



1 INTRODUCTION 

The Fe Ka spectral complex is the most prominent spec- 
tral feature in a variety of sources observed with X-ray ob- 
servatories such as ASCA, Chandra, and XMM-Newton . 
Whereas the singular 6.4 keV feature is thought to arise 
from a relatively 'cold' plasma through K-shell ionization 
by hard X-rays and subsequent fluorescence in closed 2p- 
subshell systems (Fe I - XVI), the higher energy features at 
6.7 and 6.9 keV would arise from 'hot', highly ionized He- 
like Fe XXV and H-like Fe XXVI. Other ionization species 
of Fe may also contribute in between these energy ranges. 
It is well known that in the highly ionized regime the He- 
like spectrum often dominates X-ray spectra since the closed 
K-shell survives in the plasma over the largest tempera- 
ture range and as such is a powerful diagnostic of tem- 
perature, density, and io nization state (^Gabne^^JoM^ 
19691: iMewe fc Schriiveil IT978a.b: Pradhan fc Shulll Jl98l|: 
Pradhanll98lll985bl:|Liedahi200(]l:IPorauet fc Dubaul20od: 
Bautista fc KallmanI [20001 : iNess et alJ |2002^ . The He-like 
spectra are of equal interest in laboratory plasmas such 
as tokamaks, magnetic Z-pinch and inertial fusion devices 



where non-equilibrium, time-dependen t ^ and non-thermal 
cond itions prevail (e.g. iLee et all Il985l . Il986l : iBailev et alJ 
|2002). 

In astrophysical situations the variation of the contin- 
uum and the line flux reveals infor mation on the intrin- 
sic nature of the source. For example, lYaooob et al.l ll2003l) 
have reported time-resolved Chandra and Rossi Timing X- 
ray Explorer (RXTE) detection of Fe XXV and Fe XXVI 
Ka complexes from the narrow-line Seyfert 1 galaxy NGC 
7314 and found that the ionization states of Fe respond to 
continuum variations. While the ambient plasma may still 
be in quasi-equilibrium, the flare model for the irradiated 
accretion disc is invoked to explain the Fe K a feature and 
related variations in X -ray and UV fluxes dPounds et alJ 
I2OO3I: ICollin et al.ll20o3) . Another example is that of time- 
dependent X-ray aft erglows i n gamma-ray bursts showing 
Fe Ka features (e.g iRees fc Meszaroa .200(tl . A flare or a 
burst might ionize through a number of ionization states, 
and since recombination time-scales are much longer than 
the duration of the event itself the spectrum should reveal 
the state of ionization. The models must therefore consider 
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all ionization states, particularly those likely to dominate 
under most conditions. 

Several recent observations of variable spectra of the 
iron Ka complex in AGN sho w that it is primarily due to 
Fe XXV ilwasawa et al.ll2004) . The XMM-Newton observa- 
tions of the prototypical narrow-line Seyfert 1 galaxy I Zw 
1 show an X-ray flare and spectral v ariability with both the 
Fe X XV Ka feature and the K edge llDewangan fc Griffiths! 
|20o3)- While currently observed Fe XXV spectra are unre- 
solved, and unable to constrain the associated parameters 
precisely, a more detailed analysis based on higher resolu- 
tion measurements than possible from Chandra or XMM- 
Newton might be possible in the near future. Such a pos- 
sibility should be afforded by the upcoming mission Astro- 
E2 which aims to resolve the Ka complex of Fe XXV. The 
present work aims to facilitate such studies. 

More generally, astrophysical plasmas with He-likc iron 
vary widely. In the solar corona one has the characteris- 
tic 'coronal' or coUisional equilibrium where photoioniza- 
tion may be neglected. The standard coronal spectrum of 
Fe XXV in equilibrium at different temperatures was re- 
ported in an earlier work, including the principal and di- 
electroni c satellite lines and relative contributions to the Ka 
complex dOelgoetz fc PradhanllioOlh . The present work is a 
more extensive study of the time/density dependence of a 
high temperature plasma as reflected in spectral variations 
within the 6.7 keV complex of Fe XXV. While all ionization 
states of Fe are included in the time-dependent simulations, 
we concentrate on the detailed features of Fe XXV since 
it is likely to be the prime contributor to Fe Ka in any 
highly ionized plasma. The results are obtained using a new 
code HELINE designed to study He-li ke ions in general em- 
ployi ng state-of-the-art atomic data iOelgoetz fc PradhanI 
120011) . A comparison of the non-equilibrium spectra is made 
with both coUisional and/or photoionization equilibrium as- 
suming different electron densities and radiation fluxes. Ion- 
ization fractions and corresponding spectra are computed 
with time-temperature profiles related to observed sources. 



2 THEORY 

The theory of spectral formation by H e- like ions has been 
develo pe d bv [ Gabriel fc JordanI ^1969^: [B lument hal et alJ 
J1972D: ITinetal. ! lll97iD: iMewe fc Schriiyeij <ll^78a^ : 



has al so been discussed rece n tly, am o ng others, by|Liedah 
i 2000[) : IPorauet fc Dubaul i2000h : iBautista fc Kallmar 
1 20001) . 



iPradhan fc ShulJ Jl98lD : Il3elv-Dubau et alJ JlQsl) . These 
studies assumed stationary ionization equilibrium condi- 
tions. Emissio n from non-stationa r y trans i ent p lasmas was 
considered by iMewe fc Schriiveil lll978bl . Il980l) : however, 
atomic data employed in their work has been superseded 
by more recent calculations (described later). In particular, 
the excitation rate coefficients used in their work were com- 
puted without the important effect of resonances on cross 
sections. Employing improved excitation rate coefficients 
for He-like i ons including resonanc es in the close cou pling 
formulation JPradhan et al.l | l981allbt |Pradhanlll983allbl) . it 
was shown bv lPradhanI (|l982,ll98S3) that, in particular, the 
triplet/singlet ratio G (see the definition below) depends 
considerably on the temperature and the ionization state 
of the plasma. Laboratory stu dies of ionization and recom- 
binat i on do minated plasmas (iKallne et al ] |l984l: iLee et alJ 
Il985l. ligSeT) verified the expected behavior of G(T). This 



2.1 He-like iron 

Earlier works on He-like iron have desc ribed different as- 

J ects o f spectral variations. For example. iMewe fc Schriiveil 
1978d) considered the time dependence of line ratios and 
ionization fractions, but not the d etailed spectrum includ - 
ing DR satellites. In a later study iMewe fc Schriiveilll980h 
they examined time dependent spectra of a number of ion- 
ization states of Fe and Ca under solar flare conditions, 
without an external radiation fie ld. The effect of photoex - 
citation was first p o inted out bv iBlumenthal et alj il972l) . 
ISwartz fc SulkanenI 1^9^ calculated the full Ka complex, 
but not the time dependent behavior of individual features. 
The present numerical simulations include all the features 
mentioned using new atomic data. 

Recently, emission from the 6.7 keV complex of He- 
like Fe XXV has been analyzed in ionization equilib- 
rium including the important contribution by dielectronic 
satellites to the Ka complex at T < 6.7 keV possibly 
preva lent in photoionized plasmas iOelgoetz fc PradhanI 
I2OOII) . They consider the set of forbidden (f), intercom- 
bination (i), and resonance (r) lines, with the proper 
spectroscopic designation labeled as w,x,y,z corresponding 
to the 4 transitions to the ground level Is^ C ^0) < — 
ls2p(^Pi°), ls2p(^P2°), ls2p(^P°), ls2s(^5'i) respectively. In 
recombination dominated, non-coronal plasmas (e -I- ion) re- 
combination preferentially enhances the triplet (x,y,z) lines, 
and hence the ratio G = (x-|-y-|-z)/w depends on the effective 
level-specific recombination rates, and therefore the H-like to 
He- like ratio (Xpg xXVf/^Fe XXv)' Generally, G « 1 ^ 
coUisional (coronal) ionization equilibrium; G > 1 ^ recom- 
bination dominated (photoionization or hybrid); G < 1 ^ 
ionization dominated. 

But in addition to the principal lines the dielectronic 
satellites (DES) are observe d to be quite prominent in the 
laboratory ( e.g. iBei crsdorf cr et al.lll992^ . solar flares and 
corona (e.g. iKato et alJll99d) . Therefore the Ka complex 
of He-like Fe XXV needs to be studied using a model 
that includes the n = 2 bound levels, as well as the 
KLn autoionizing levels of He-like ions. Using our ear- 
lier time-independent model in the HELINE code. Fig. 1 
from lOelgoetz fc Pradheinl (|20^) shows the Ka complex of 
Fe XXV computed in coronal equilibrium at several temper- 
atures. At temperatures well below th at of maximum abun- 
danc e of Fe XXV, T < T^ « lO'^ * K jArnaud fc Ravmondl 
Il992h . the DES dominate Ka emission. On the other hand, 
for T >Tm the DES intensities are small relative to that 
of principal lines w,x,y,z, with large variations in between. 
In lower temperature plasmas (e.g. j<^allnc ct a]^ (.1984}) ) and 
photoionized sources we expect the centroid of the Ka com- 
plex of Fe XXV to lie shortward of the w-line and show 
'red ward broadening' of the ot herwis e 6.7 keV feature. 



Swartz fc Sulkane: 



also 



m coro- 



Jl993f) (see 
lOclgoctz fc PradhanI I2OO j) further show that 
nal equilibrium the ratio G(T) for Fe XXV is approximately 
a constant, independent of T. However, it is clear from Fig. 
1 that the n = 2 KLL DES would be largely blended in 
with the principal lines since many of the strongest ones 
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lie in between the z and the w hnes. Therefore, foUowin R 
several earher works su c h as [ Swartz fc SulkanenI il993l) . 
iBelv-Dubau et al.l ^^E^; iKato et al.l l|l99Sh . Oelgoetz and 
Pradhan assume that the DES(KLL) intensity is merged 
with the (x,y,z) lines (see Fig. |5| ) , but may be measured 
relative to the w-line blended with the higher-n > 2 
satellites. They they re-define the G ratio, including DES, 
as 



GD 



{x + y + z + KLL) / (w + KLn) 



n>2 



(1) 



lOelgoetz fc PradhanI (1200 if) established the tempera- 
ture variations of the Ka complex in equilibrium. They also 
plot GD(T) from equation Q as a function of T and the 
ionization fraction Xi/ . It is ascertained that GD(T) falls 
sharply at T >Tm,, owing to the decrease DES emission. 
Nonetheless, the effect of departure from ionization equi- 
librium was estimated in the time-independent model in 
an ad hoc manner. In this study we introduce temporal- 
temperature dependence explicitly. 



2.2 Time-dependent plasmas 

It is generally assumed that ionization/recombination rates 
are much slower than radiative-coUisional rates. Therefore 
the temporal-temperature (t,Te) variations would manifest 
themselves via the ionization fractions Xi obtained from the 
coupled ionization balance equations 



dXi 
dt 



+ar+'l)X,+i + C.-i(T, 



hf 



Gi{hi')dv 



+Cf"{T) + Cf^{T) 



(2) 



(3) 



where af ^ and af^' respectively are the dielectronic and 
radiative recombination rates out of ionization state i, C; 
are the rate coefficients for collisional ionization, with Cf ^ 
and Cf^ as the direct and excitation-autoionization con- 
tributions respectively, (f){hv) the radiation energy density, 
and cTi the photoionization cross section of ionization state i. 
Three body recombination is neglected, as it is unimportant 
for densities presented in this study. Note that the ionization 
fractions would vary with the transient behavior of any of 
the three parameters: Te, Ne, or (f>(hu). The radiation field 
intensity may be expressed in terms of the conventional ion- 
ization pa rameter U, th e ratio of the local photon to electron 
densities jKroliklll999ll 



U 



1 



l.ORyd 



hv 



dv 



(4) 



Once the ionization balance equations are solved, and 
the fractional populations of each ionization state found, 
the level populations of the He-like ion are then obtained by 
solving 



-7VeiV.^gr^+iVe5]iV,g;:(r) 
-N,Ni J2 1?j + NeCt'{T)Xu 
-N,Ch{T)N, + N,(j>{hPji)Bji 
-N,Y(t){hu,j)Bij 



(5) 



where the A's are the radiative decay rate coefficients, B's 
absorption and stimulated emission rate coefficients, Qi 's the 
level specific recombination rate coefficients, C's the colli- 
sional ionization rate coefficients, qji, q"i, and qfj the rate 
coefficients for excitation/dexcitation by electrons, protons, 
and a particles respectively. 

The problem is then further simplified by using effec- 
tive rates, including cascade corrections to redu ce the model 
to a seven level problem. Unified level specific (iNahar et alJ 
2001^ rates, are supplemented with distorted wave type ra- 
diative recombination and dielectronic recombination rates 
calculated using the Los Ala mos edition of Cowan's Atomic 
struc ture code, CATS (see: ICowanI Il98ll : lAbdaUah et. al 



W8i ) and the ionization code GIPPER (seeT f^reh^^t^^ 



20001) where the more detailed unified level specific rates 



are more difficult to obtain accurately. Electron impact ex- 
citation rates into and out of high lying states (states for 
which close coupling calculations including resonances are 
unavailable) are also calculated in a similar manner, again 
via the distorted wave app roximation, using the Los Alamos 
collisional code ACE (see: IClark et al.lll988f) . 

In addition, the satellite line intensities due to dielec- 
tronic recombination (DR) and inner shell excitation (from 
the ground st ate) are calculated using the expressions of 



tDR . 3/2 3 V- ^^ Tn-3/2 -1-^ QaVrVa 



Ta + E 



Ta + E 



(6) 
(7) 



where Tr is the radiative decay rate coefficient. Fa the total 
autoionization rate coefficient, Es is the energy of the tran- 
sition, Qa the statistical weight of the upper state, and q{Te) 
is the electron impact excitation rate out of the ground state 
and into the autoionizing state that produces line s. 

Inner shell photoexcitation is taken care of in a straight 
forward manner using the expression: 



= XLi4>{hVs) 



Fa + E F, 



(8) 



using the same variable definitions above. 

It should be noted that all coupling between the au- 
toionizing states that give rise to these satellite lines is ne- 
glected. 
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3 COMPUTATIONS 



As discussed above, first we solve the coupled ioniza- 
tion balance equations and |^. CoUisional ionization 
along with both dielectronic and radiative recombina- 
tio n rates are calculate d on t he fly using the expressions 
of lArnaud fc Ravmondl il992l) . Photoionization rates are 
pretabulated for a number of radiation fields by integrating 
over t he resonance averaged cross sections of iBautista et all 
il998h and then read in. If a steady state solution is de- 
sired, the matrix representation of the coupled ionization 
balance equations is factored via Gaussian elimination and 
then solved by the standard LINPACK routines. Time de- 
pendent ionization fractions are found by employing the 
Fehlberg fourth-fifth order Ilunge-K utta ijFehlberd Il96fll ') 
routine bv lShampine &i Wattj (Il977fl . The initial ionization 
fractions are either provided to HELINE in the input file, or 
more commonly the initial ionization fractions are those of a 
steady state plasma, and as such are calculated. The result- 
ing equations are then integrated forward in time, allowing 
the routine by Shampine and Watts to handle the choice of 
time step. 

Once the ionization fractions have been calculated for 
a given time or set of steady state conditions, spectra can 
be produced by solving the coupled set of equations that 
arise from equation|K| We assume that the time scale for the 
ionization and recombination rates is much longer than the 
time scale for the radiative rates, as such equationl^is always 
solved in a ti me independent manner . Th e collisional excita - 
tion rates of IWhiteford et al.1 (1200 ill and lPradh aTi figSSaf), 
the lev el specific unified recombin ation rates of j Naha r ct alJ 
i200lfl and radiative dec ay rates of[Nahar fc Pra^rMTin^ggf ) 
are used along with the lLota formula for inner shell 

ionization. Furthermore, ionization and recombination out 
of the n=2 level s of the He-like species is neglected, follow- 
ing the model of lMewe fc Schriiveil lll978af) . Again Gaussian 
elimination is employed to first factor the matrix representa- 
tion of the system of coupled equations, and then the prob- 
lem is solved using the LINPACK routines. The population 
of each of the four states that give rise to the four He-like 
lines of interest (w, x, y, z) are then multiplied by the cor- 
responding radiative decay rates to produce line intensities. 
Satellite line intensities are calculated using equations 1^1 171 
|S| T he autoionization a n d rad iative decay rates are taken 
from IPradhan fc Zhand Jl997t) : the inner shell excitation 
rates of nCat^^^dTi^gsl) are used as well. 

Once line intensities have been calculated, each line is 
given a Gaussian line shape (with a width determined by 
the electron temperature). Component spectra arising from 
radiative decay of the He-like spectroscopic states, both in- 
ner shell electron impact excitation and photoexcitation of 
the Li-like ground state, and dielectronic recombination into 
the KLL autoionizing states of the Li-like ion are written out 
separately, along with their sum, the total spectra. The se- 
quence then repeats, either with the next steady state ion- 
ization balance calculation, or by integrating the time de- 
pendent ionization balance equations forward to the next 
point of interest. 



4 RESULTS AND DISCUSSION 

We compare and discuss steady state and transient plasmas 
with respect to (i) ionization balance, and (ii) the 6.7 keV 
Ka spectral features. Their relevance to the standard He 
quartet line ratios R = z/(x-|-y) and G = (x-fy-|-z)/w are 
pointed out (It is usual practice to denote the forbidden (f 
or z), intercombination (i or x-l-y), and resonance (r or w) 
lines of He- like ions as a 'triplet'; however, for highly charged 
ions such as Fe XXV where the fine structure separation be- 
tween x and y is manifest and the excitation atomic physics 
is distinct for each line, it is more appropriate to refer to 
them as a 'quartet', as herein), but more general diagnos- 
tics are developed including (a) the effect of radiation fields 
on ionization balance and photoexcitation of levels, (b) high- 
density effects, and (c) spectral time- variations in ionization, 
equilibrium, or recombination dominated plasmas. 



4.1 Steady State Plasmas 

Under collisional or photoionization equilibrium the ioniza- 
tion balance may be fixed as function of temperature. It is 
then straightforward to obtain the spectral intensities using 
HELINE with dXi/dt = in Eq. (3). 



4-1.1 Ionization Balance 

As described in the Computations section, the first step in 
any calculation is to solve the coupled set of equations aris- 
ing from equations|5|and|3 Results of this step as a function 
of electron temperature for two different radiation fields, a 
blackbody and a power-law, along with a purely collisional 
case (J = 0), are presented in Fig. The radiation fields 
chosen are purely exemplary, after a number of trials, so as 
to highlight typical features. Each radiation field is deter- 
mined by either a Planck function or a power-law with a 
specified photon index. 

The first thing one might note is that the radiation field 
allows a given ionization state to become abundant at lower 
electron temperatures relative to the pure collisional case. 
Additionally, the range of electron temperatures over which 
a given ionization fraction can be found can be significantly 
broadened. It should also be noted that for a faster decay- 
ing radiation field, such as a blackbody, the radiation field 
tends to effect the lower lying ionization states, but not the 
most highly charged ionization states that depend most sig- 
nificantly on the high energy tail. This effect can be seen by 
comparing the ionization fraction curves for Fe XXIII - Fe 
XXVII in the upper (blackbody) and the lowest (pure col- 
lisional) panels of Fig. Q with the middle panel (power- law 
with a photon index a — -2.3). Lastly it should be noted 
that by working in terms of U (dividing the radiation field 
by A^e) all the ionization and recombination terms are first 
order with respect to electron density (A'^e) and as such the 
ionization balance curves (as a function of electron temper- 
ature) in equilibrium only depend directly on the value of 
U and the shape of the radiation field, not on the over all 
intensity of the radiation field, or on the electron density. 

The present approach, however, is explicitly predicated 
on the assumption that in some transient plasmas ioniza- 
tion equilibrium or energy balance is not satisfied globally 
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or locally. For example, fusion plasmas with sudden exter- 
nal energy input, or in flaring situations in astrophysics, may 
not be constrained by parameters such as direct coupling of 
temperature and radiation fields, cases normally dealt with 
conventional photoionization equilibrium codes. This is one 
of the crucial differences of the present method with exist- 
ing codes; the present code HELINE is therefore applicable 
to transient phenomena in both the laboratory and astro- 
physical situations. The value of the ionization parameter U 
is computed a posteori to illustrate correspondence with a 
local radiation and electron density, but is not confined to 
them. In fact, no single temperature may be defined that 
corresponds to an exact particle or radiation field distribu- 
tion. HELINE attempts to simulate this decoupling in a va- 
riety of ways, while also representing photoionization and/or 
collisional equilibrium. 

4.1.2 Spectra 

The first patterns shown in the spectra of Fig. |5| (panels in 
the first 4 rows) are the approximately linear scaling of in- 
tensity with electron density, and the importance of satellite 
lines at lower temperatures (LHS column) . The former is due 
to the fact that collisional excitation, collisional ionization, 
and recombination are the principle mechanisms of produc- 
ing ions in the excited state, and all are first order with 
respect to electron density. The intensities of Li-like satel- 
lites produced by inner-shell collisional excitation is greater 
at lower temperatures for two reasons. First one can see in 
equation |n] that electron temperature appears both in the 
denominator, and in an exponential term. As the electron 
temperature decreases, the dielectronic satellite lines grow in 
strength owing to the dominance of the exponential term in 
this temperature range. Secondly, the Li-like ionization frac- 
tion appears in equation |7| and Li-like iron ionizes quickly 
as one can see from Fig.Q 

In order to simulate the effect of a relatively intense ra- 
diation field, such as may correspond to flaring conditions, 
we choose a fairly large value of U ~ 200. The results are 
shown in the bottom two panels of Fig. Q . The basic effect 
is analogous to that due to high densities (discussed in the 
next section), i.e. to redistribute level populations. However, 
the detailed effects of radiative and collisional transitions are 
different. Radiative transitions strongly favor dipole allowed 
ex citations with large oscilla tor strengths. It was first noted 
by iBlumenthal et al.l that the 2^5 2^P transi- 

tion is the most likely to be affected by photoexcitation. It 
should be noted however that these dipole allowed absorp- 
tion, and stimulated emission transitions must compete in 
overall rate with the radiative decay rates in order to be 
important, thus unlike the ionization balance problem, U is 
not the important parameter. In the case of the transitions 
mentioned above, this implies strong radiation fields for all 
electron densities. 

In order to keep the U values, and hence the ionization 
balance in a physically reasonable state, the increment in 
the intensity of the radiation field was done in conjunction 
with an increase in electron density. Of the four transitions 
mentioned above, the transition from the 2^ Si level to the 
2^P2 level is preferred as it has the greatest B-coefficient. 
The effects of photoexcitation along the 2"^ Si 2^Po tran- 
sition would not ordinarily be visible, as there is no spec- 



tral line produced by the decay of the 2^Po level. The other 
transitions have similar B-coefflcients to each other, and as 
such the transitions 2^Si 2^Pi,2 are radiatively excited. 
As one would expect, this results in a diminishing of the 
z line, and an enhancement of the intensities of predomi- 
nately the X and to a lesser degree the y line. The net effect 
is that photoexcitation decreases the usual electron-density 
sensitive line ratio R = z/(x-|-y), while G is essentially un- 
changed since it is the ratio of the sum of all triplets to the 
singlet w line. 

Figure|5|also shows the effects of electron density on the 
spectra. The general trend is that as A''e increases, popula- 
tion is moved from the metastable 2"^ Si level to the higher 
levels within the n = 2 complex. As the 2^ Pi, 2 levels are 
closer in energy than the 2^ Pi, the first observed effect is 
the transfer of population into the levels and, as its de- 
cay rate is faster, a preferential enhancement of the y line 
over the x line. While at lower densities the x line increases 
faster than the y line owing to the higher collisional excita- 
tion rate of 2^P2 than 2^ Pi from the 2^ Si, as the electron 
density increases the n = 2 complex moves closer to statis- 
tical equilibrium. Since the y line (2^ Pi l^So) and the w 
line (2^Pi — > l^So) have significantly faster decay rates than 
the X (2''P2 — > l^So) line, the x line will weaken relative to 
the others as equilibrium inside the n = 2 complex is ap- 
proached. These collisional processes at high densities even- 
tually render the z line undetectable as can be seen in Fig.|21 
Again the net effect is a decrease in the value of the line ratio 
R while leaving G essentially unchanged; also the behavior 
is different from the photoexcitation case. Should the elec- 
tron density continue to be increased, the levels would finally 
come into thermodynamic equilibrium which would quench 
the X line as well. In order to calculate accurate spectra at 
these densities however, our model would need to modified to 
include three-body recombination, and other plasma effects 
such as line broadening, which would become important for 
higher densities {Ne > 10^"). Nevertheless, the main effect 
of z-line quenching at high densities is amply demonstrated 
in Fig. Inland is a potential plasma diagnostic of sources such 
as laser produced fusion plasmas. 

At sufficiently high densities the DES will also be af- 
fected, as collisional rates among closed spaced autoioniza- 
tion levels approach typical autoionization rates of 10^''"^'' 
sec~^ at electron densities in excess of 10^° cm~^. We do not 
consider such regimes herein, but further work is in progress. 

4.2 Transient Plasmas 

Characterizing the spectral 'signatures' of transient or non- 
steady state plasmas is the primary goal of this study. We 
assume sample asymmetric time dependent functions for 
electron temperature and radiation field intensity with a 
faster rise and a slower decay (such as under fiaring con- 
ditions), shown in the upper two panels of Fig. |21 For il- 
lustrative purposes we assume a constant electron density 
(A'^e) of 10*cm~^, which is appropriate for the chosen time 
scale of a few thousand seconds (a time scale such as what 
was observed bv .Dewangan fc Griflith£..200 j) . If the density 
is much lower, the plasma is unable to respond on the time 
scale of a few thousand seconds. If the density is much higher 
the plasma will tend toward equilibrium much more quickly. 
Not only does this make for a quasi steady state plasma 
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(which was not our objective) but it also becomes incredi- 
bly time consuming to solve in a time dependent manner as 
the time steps needed to accurately solve the ionization bal- 
ance equations get inordinately small. Three time dependent 
models were used. In the first just electron temperature was 
varied (as in the top panel of Fig.|^; there was no radiation 
field. This is the coUisional case; it's ionization fractions can 
be found in the middle panel of Fig.j^and it's spectra in the 
left column of Fig. 0] The radiation field intensity was varied 
in the second case (as in the second panel of Fig. |1J and the 
electron temperature was held constant at Te = 10^ °. This 
is the photoionized case; it's ionization fractions are in the 
fourth panel of Fig.|3 and it's spectra in the middle column 
of Fig. m Both elect ron temperature, and radiation field in- 
tensity were varied in the third case. This hybrid case give 
rise the the ionization fractions found in the bottom panel 
of Fig. 13 and it's spectra in the right most column of Fig.|l| 
The results and diagnostics derived from these three cases 
is discussed in the next three subsections. 

^,2.1 The Transient CoUisional Plasma 

Ionization balance in the transient plasma is dominated by 
ionization and recombination at different times, as shown 
in the middle panel (2nd from Top) of Fig. El and the left 
column of Fig. 21 Initially (about the first 1500 seconds), 
the inner shell coUisional excitation (ISCE) lines are very 
prominent due to an overabundance of Li-like Fe XXIV and 
a higher than usual electron temperature - hence higher col- 
lisional excitation rates. The z line is also enhanced due to 
inner-shell coUisional ionization (ls'^2s + e^ — > ls2s(^Si)-f 
2e~) during the early part of this ionizing regime. The ef- 
fect on the line ratios is a larger value of G(GD) than one 
would otherwise expect for a given temperature, due to en- 
hanced z and satellite lines. In addition the G(GD) line ra- 
tios would be large because the w line may not be excited 
very early during the ionizing phase. The ratio R would be 
similarly affected, although resolving the x and y lines might 
be difficult in this case. It might be noted that this scenario 
refers to astrophysical or laboratory plasma sources where 
the coUisional and radiative rates are averaged over distri- 
bution functions: Maxwellian for particles and Planckian for 
the radiation field. The resultant spectra may be different 
in controlled experimental setups where the incident elec- 
tron beam can be mono-energetic and tuned to scan through 
the n = 2 levels. Selective excitation b y such a beam (as 
in El ectron-Beam-I on- Tra ps (EBIT), e.g. iBeiersdorfer et alJ 
Il992l: iTakacs ct aQ 120031) may show a different signature 
during the ionization phase via tuned inner-shell ionization 
of the Li-like ions. The EBIT beam widths may be repre- 
sented by a Gaussian whose widths are sufficiently mono- 
energetic to excite a given set of levels primarily by elec- 
tron impact excitation. Such simulations have been done 
for line r atios measured in EBIT spectra of Ne-like iron 
Fe X VII JChen fc Pradhanll2002l : fchen. Pradhan fc EissneJ 
l2003f) . using new Breit-Pauli R-matrix (BPRM) coUisional 
cross sections, and rates computed with several electron dis- 
tribution functions. Similar computations are under way for 
line ratios for several He-like ions using new coUisional and 
radiative BPRM data. As such, the EBIT experiments pro- 
vide another means of line diagnostics in plasmas different 
from astrophysical sources. 



In middle panels in the LHS column of Fig. 2] the tran- 
sitional period produces spectra that show a combination of 
recombination and ionization, and look like the steady state 
spectra of Fig. |5| Indeed, at ~ 1560s, the spectra reflect 
quasi steady state conditions with similar line intensities as 
in Fig. 121 although some differences such as enhanced inner 
shell excitation lines due to an overabundance of the Li-like 
ionization state are apparent. 

As recombination time scales in transient sources are 
likely to be longer than ionizing time scales, the most com- 
mon case occurs in the late phases (bottom panels of Fig.^J 
when one can see the switch to a regime where recombina- 
tion begins to dominate the plasma conditions. Several fea- 
tures manifest themselves, driven mainly by recombination- 
cascades to the upper levels of the n = 2 Ka complex. (I) 
The non-dipole x line becomes stronger than the y inter- 
combination line (>1600 seconds). (II) The forbidden z line 
is also enhanced in intensity. (Ill) As one can see from Fig. 
Olthe pi asma becomes more and more ionized relative to the 
steady state as the electron temperature continues to cool 
faster than the plasma can recombine. This lagging forces 
the plasma into a regime where recombination will domi- 
nate spectral formation. (IV) As recombination dominates 
over coUisional ionization and excitation, the level popula- 
tions in the P2 and the "^Si levels are increased relative 
to others since recombination-cascades preferentially popu- 
late the triplet levels jPradharn . 1 982,. , 1 985bll . This resuUs 
in the enhancement of the z and x lines that can be clearly 
seen from about 1900 seconds on wards. The effect is so pro- 
nounced that the forbidden lines z and x actually become 
stronger than the two dipole allowed lines w and y. We then 
have a 'spectral inversion' compared to the steady state or 
ionizing conditions shown in the other panels of Fig.|l] This 
inversion would lead to an enhanced G(GD) ratio, and to 
an R value less than 1. 

4-2.2 The Transient Photoionized Plasma 

The fourth panel of Fig. Inland the middle column of Fig. 21 
shows results for a sample, transient, photoionized plasma. 
There are similarities with the coUisional model in that it 
has distinctive ionizing and recombining phases and that 
the transitional stage between them blends the features, but 
the most outstanding aspect of photoionized plasmas is that 
throughout the process the DR and ISCE satellite lines are 
weak; the exponential terms in eq. |H1 and |7| dominate and 
are too small to allow the satellites to form. Thus only the 
narrow satellite lines from photoexcitation (q, r and t) along 
with a weak y line and a strong w line produce the spectra in 
the ionizing phase. This is a huge difference from the afore 
described coUisional case (in the panels of the left column). 
This results in a blue-ward shift of the average flux and a 
G(GD) ratio which is much closer to 1 for much of the early 
ionizing phase. 

The plasma then begins a gradual transition, as recom- 
bination becomes more important with increasing time. At 
~1080 seconds the plasma is just entering the recombination 
phase, as is evident by the now stronger non-dipole lines 
populated by recombination-cascades. The line intensities 
(in this case the He-like line intensities only) are compara- 
ble and the spectra looks qualitatively like equilibrium spec- 
tra, but without satellite lines. At later times however, the 
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plasma becomes dominated by recombination, and by about 
1900 seconds the z and x lines have become stronger than 
the dipole allowed w and y lines. By about 2400 seconds the 
spectra of this transient photoionized case looks very similar 
in shape to the spectra of the transient coUisional case. 



4.2.3 The Transient Hybrid Plasma 

When both the electron temperature and the radiation field 
intensity are varied in a time dependent manner, one would 
expect to see a mixture of the features of the previous two 
cases. The bottom panel of Fig.j^and the right hand column 
of Fig. 2] shows such a situation. 

First the plasma passes through an ionizing phase t < 
1500 seconds. During this time the w line is markedly en- 
hanced along with the satellite lines produced by inner shell 
photoexcitation; the z line, along with the satellite lines 
arising from inner shell coUisional excitation, are stronger 
than in the photoionized case, but not dominate as in the 
purely coUisional case. This indicates that both coUisional 
and photoionization are ongoing, although photoionization 
is the more important of the two, as can be seen by the 
prominence of the w and q lines. One need only look at the 
ionization fractions presented in the bottom panel of Fig. |H] 
to see that the radiation field presented is more effective at 
ionizing than the ramping electron temperature used. 

By t ~ 1560 seconds the plasma has made the transition 
into a recombining mode, although it has not yet begun 
to show the characteristic spectral inversion. One note of 
interest is that because the electron temperature is higher 
than in the transient photoionized case, the DES lines j and 
k can be seen, even if the principal He-like lines are more 
intense. By t ~ 1920 seconds the spectral inversion that is 
common to all three models at late times is evident. 



5 CONCLUSION 

High resolution steady state and time dependent analysis of 
the Ka complex of the primary (w, x, y, & z) lines and di- 
electronic satellites is carried out in unprecedented detail. It 
spans an energy range between 6.6-6.7 keV; with a signature 
sensitive to ionization, recombination, temperature, density, 
and radiation field. 

The salient features of the simulated spectra are: 

(i) The dielectronic satellites are prominent around 6.65 
keV at low temperatures, such as in some photoionized plas- 
mas, shifting the Ka complex intensity red- ward to energies 
below the w line at 6.7 keV. In general lower energy features 
may be present during both the ionizing and recombining 
phases, although the latter is more likely. 

(ii) Photoexcitation and electron density sensitive level 
redistribution affects the x and y lines differently; pho- 
toexcitation tends toward a stronger x line, while electron 
density towards a stronger y line, and as such they may 
be able to be distinguished if the spectra are adequately 
resolved. 

(iii) A time dependent source such as an X-ray flare can 
induce a number of observable spectral changes: 



(a) Enhancement of the z line and select satellite lines 
(q, r, t, u and v) due to inner shell coUisional excitation 
lines during coUisional ionization, resulting in higher than 
usual G(GD) ratio in particular. 

(b) Enhancement of the resonance (w) and select satel- 
lite lines (primarily q, r, and t) during photoioniza- 
tion, resulting in a lower G(GD) ratio before the H-like 
Fe XXVI ionization fraction, and recombination-cascades 
therefrom, becomes important. 

(c) Enhancement of the non-dipole (z and x) lines and 
the 'spectral inversion' with respect to the dipole allowed 
y and w lines during recombination. 

(iv) In order to realize the full potential of the line diag- 
nostics developed in this work, high resolution observations 
are needed. This should be possible in X-ray spectroscopy 
from Astro-E2 that is expected to resolve th e Fe XXV 6.7 
keV complex with up to 6 eV resolution iYagoob et alJ 
l2003h . The basic diagnostics should also be applicable 
to inertial confinement fusion devices with nanosecond 
time scales, but considerably higher densities than in 
astrophysical sources. Such simulations are in progress. 
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Figure 1. Steady state ionization fractions for three different radiation fields as a function of electron temperature: top panel a blackbody 
{Bhv(Thu = lO^'^)), middle panel — power-law with a = -2.3, bottom panel — coUisional equilibrium without radiation field. Both 
radiation fields correspond to same ionization parameter U ~ 200. Notice that while the power law is not as effective at photoionizing 
lower ionization states, it has a greater effect on the higher ionization states due to it's slower decay. 
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Figure 2. Steady state spectra at a variety of densities and radiation fields. In particular, the spectra demonstrate the quenching of the 
forbidden z line under different conditions. All intensities are given in the same set of arbitrary, relative units. 
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Figure 3. Time dependent profiles: Top: the electron temperature, 2nd from Top: ionization parameter U(t) associated with the radiation 
field, 3rd from Top: ionization fractions under coUisional ionization alone (Te(i) as above), 2nd from Bottom: ionization fractions for the 
transient photoionizod case (variation of the radiation field intensity, U(t) as above and constant Te = 10®"), and Bottom: ionization 
fractions from the hybrid case (varying both Te(t) and U(t)) 
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Figure 4. Evolution of transient spectra at various times throughout the time-dependent runs, with Te{i) and U(t) as in Fig. |^ (relative 
intensities). The LHS panels correspond to pure coUisional case, the middle panels to photoionized case, and the RHS panels to the hybrid 
case. Of particular note are the enhanced satellite and z lines at early times in the coUisional case (signature of coUisional ionization), the 
enhanced resonance (w) and inner shell photoexcitation (q,r & t) lines at early times in the photoionized case and hybrid cases (signature 
of ongoing photoionization), and the enhanced non-dipole (x and z) lines at late times in all cases (signature of recombination). All 
intensities are given in the same set of arbitrary, relative units. 
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